Many-body Rabi oscillations of Rydberg excitation in small mesoscopic 

samples 
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We investigate the collective aspects of Rydberg excitation in ultracold mesoscopic systems. 
Strong interactions between Rydberg atoms influence the excitation process and impose correla- 
tions between excited atoms. The manifestations of the collective behavior of Rydberg excitation 
are the many-body Rabi oscillations, spatial correlations between atoms as well as the fluctuations 
of the number of excited atoms. We study these phenomena in detail by numerically solving the 
many-body Schredinger equation. 

PACS numbers: 32.80.Rm, 03.67.Lx, 34.20.Cf 



Introduction 

Long radiative lifetimes of Rydberg atoms and 
their possibility to interact strongly at large dis- 
tances have made ultracold Rydberg atoms in- 
teresting systems for possible quantum informa- 
tion applications. The electron in a Rydberg 
state is very far from the nucleus and thus sen- 
sitive to external fields or the presence of neigh- 
boring Rydberg atoms. Due to the huge polar- 
izabilities of Rydberg atoms, it is possible to in- 
duce relatively large electric dipole moments us- 
ing small electric fields. The capability to turn 
on and off the interactions just by switching the 
external field is an important aspect of this ap- 
proach to quantum computation. In this way 
decoherence effects due to the interactions be- 
tween atoms or with the environment, can be 
significantly reduced. Strong interactions can 
be used to entangle neutral atoms and achieve 
fast quantum gates [U, 0] i as well as to blockade 
excitation by shifting many-atom excited states 
out of resonance. It has been proposed to use 
this blockading effect to realize scalable quan- 
tum gates [3| ■ The evidence of excitation block- 
ade has been found in several experiments with a 
narrow laser bandwidth. In Q, a local blockade 
of Rydberg state excitation in a mesoscopic sam- 
ple due to strong van der Waals (vdW) interac- 
tions has been observed using a pulse-amplified 
single-mode laser. The 5s ground-state rubid- 
ium atoms were excited by one photon UV tran- 
sitions to high np3/2 Rydberg states. Rydberg 
excitation exhibited dramatic suppression com- 
pared to the non-interacting case. A mean-field 
type model was proposed to explain these exper- 
imental results. In the model, different atoms in- 
teract differently depending on their locations. 



Different interaction energies were modeled by 
a distribution of mean-field shifts for which a 
distribution of excitation probabilities was cal- 
culated. A good agreement between the theoret- 
ical model and experimental measurements was 
found. Significant suppression of Rydberg exci- 
tation has also been observed using cw excita- 
tion This Rydberg excitation, strongly in- 
fluenced by interactions, exhibits sub-Poissonian 
atom counting statistics [1, 0| ■ The blockade ef- 
fect due to dipole-dipole interactions in an ultra- 
cold sample of Cs atoms 8] has been reported. 
Also, an interesting antiblockade effect in two- 
step excitation processes was predicted d, • 

At high principal quantum number n, the 
interactions between Rydberg atoms are quite 
strong and they inhibit the excitation of many 
surrounding atoms within a range of few /im. 
The atoms within this range are strongly cor- 
related so that a many-body treatment is, in 
general, needed. In [111 ] the many-body wave 
function was calculated by numerically solving 
the Schrodinger equation. This type of analy- 
sis can be irnportant for quantum information 
applications [I4I because precise control at the 
quantum level is essential in these applications. 
Strong interactions can affect the coherent ma- 
nipulation of a large group of atoms so it is nec- 
essary to use a many-body treatment to evaluate 
the fidelity of quantum information protocols. 

In this paper we consider mesoscopic sam- 
ples of ~10 1MB. diameter containing up to 100 
ground-state atoms which are excited by single 
photon transitions to high np Rydberg states. 
We study in detail the dynamics of such systems, 
especially the possibility of many-body Rabi os- 
cillations of Rydberg excitation. It is plausible 
to investigate these oscillations in smaller sys- 



2 



terns because one would not normally expect 
to achieve the coherent manipulation of large 
groups of atoms. The many-body approach de- 
veloped in [ll| is quite suitable for this analysis 
and we use it here. We only modify some tech- 
nical details on how to treat interactions in this 
approach. 
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Many-body effect in ultracold Rydberg 
systems and collective oscillations 

Although models may be very useful in ex- 
plaining some important overall properties of 
large strongly-interacting Rydberg systems, it 
is essential to include many-body correlations in 
the study of many-body effects. In ultracold Ry- 
dberg systems, the thermal motion of the atoms 
is greatly reduced so that in many situations it 
can be completely ignored. 

The idea of these collective oscillations can be 
explained as follows. If the interactions between 
atoms are strong enough, all many-body states 
with two or more excited atoms will be greatly 
shifted by these interactions and thus far-off 
resonance. Such systems are effectively two- 
level systems because there are only two (col- 
lective) states that are populated: the ground 
state \G) — |(7i(72 • ■ ■ 57v) , where all atoms are 
in the atomic ground state \gi), and only one 
excited state \E) = 1/y/N \gi . . .ei . . . g^) 
where any, but only one, atom can be excited. 
On atomic resonance such two-level systems are 
exactly solvable [3] and the solution, in terms 
of excitation probability Poxc , is 

PoAt) = ^sin^ (yNnF{t)) , (1) 

where n,F{t) is the pulse area and N is the 
number of atoms in the sample. Clearly, the 
collective oscillations are much faster then the 
isolated-atom Rabi frequency fl. The question 
is whether these fast oscillations can exist in sys- 
tems which are big enough that there can be 
several excited atoms, or alternatively with the 
interactions not strong enough to fully blockade 
the system. 

Because the atoms are ultracold, we ignore the 
effects of thermal motion during Rydberg exci- 
tation in this analysis. We consider the following 
many-body Hamiltonian of two-level atoms and 



where A is the frequency detuning and the in- 
teractions between Rydberg atoms are given by 
Kij. The second term in the Hamiltonian is 
the dipole operator representing the interaction 
with the optical field. The function f{t) is the 
time evolution (envelope) of the laser pulse. The 
non-trivial parts of the cr-operators are defined 
as al^p = \ai) {f3i\, where a, (3 reffer either to the 
ground state g or the excited state e. 

We can always eliminate the de- 
tuning from the previous Hamilto- 
nian by the unitary transformation 
exp (zt AEfliaL) H exp {~it AJ^fL.&l,) . In 
the remaining part of the Hamiltonian, the only 
effect of this transformation is to replace f{t) by 
f(t) exp(rfA). This transformation can simplify 
the use of other approximations applied to solve 
We assume that it is always performed and 
there is no need to consider this term explicitly. 
In fact, we set A = in our simulation. 

Since the dimensionality of the problem is 
2^, which is a huge number for numerical cal- 
culations with TV « 100, it is necessary to 
make some approximations. The idea of lo- 
cal blockade is utilized here (more explanations 
can be found in [11]). We can group atoms in 
such a way that the probability to have two or 
more excited atoms within a group is practically 
zero. Such groups of atoms are often called su- 
peratoms. Each superatom "i" is a two-level 
system described by two collective states \Gi) 
and \Ei). In the next step, the actual many- 
body Hamiltonian is represented/approximated 
in terms of these superatoms. For convenience 
we define new operators a^^Q — \Ei) {Gi\ and 
'^EE ~ In. terms of the new opera- 

tors, the interaction with the laser field can be 
expressed as follows 

N 

ifitK+ntKe) 

4=1 

-> E ^ (/w^^G + r , (3) 
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because 
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where iVga is the number of superatoms and Nj 
is the number of atoms in a superatom. The 
Rydberg interaction Vjiyd between superatoms 



and 



'J " in their excited states is 

J: N: 



(4) 



p<Ni q<Nj 



In these sums, atoms "p" and "q" belong to dif- 
ferent superatoms. The last formula gives the 
prescription for how to introduce the interac- 
tion between superatoms. It seems that the in- 
teraction between two superatoms in [ll| was 
modeled using the distance between their cen- 
ters of mass. Equation (jH) suggests that the 
superatom interaction is rather defined by (Kpg) 
instead of K{{rpg)). Using K{{rpq)) tends to un- 
derestimate the influence of interactions. The 
difference between these two ways can be sig- 
nificant if K{r) depends strongly on r, which is 
the case here. However, the effect of such dif- 
ferences in interaction energies on the results of 
simulations is hard to judge because the largest 
Kpq do not contribute to the interactions be- 
tween superatoms (explained later in the text) 
and because of the strong suppression of Ryd- 
berg excitation. This kij , as the relevant param- 
eter to describe the interactions between super- 
atoms, should be preferably used in the process 
of the superatom formation as well. These are 
the only modifications we make to the method 
Another useful simplification comes from 
the suppression of Rydberg excitation. Namely, 
the number of excited atoms is limited due to 
the blockade effect and there is no need to con- 
sider many-body states with high numbers of 
excited atoms. To illustrate the size of the prob- 
lem after these steps, we note that for systems 
with twenty-three superatoms and at most seven 
of them excited, there are about four hundred 
thousand excitation amplitudes to solve. 

To form superatoms we use the same recursion 
described in except that the largest kij is 
our guidance in deciding which (super) atoms to 
group together. We use the following recursion. 
We initially set the number of superatoms to be 
equal to the number of atoms. After that, we 
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FIG. 1: The number of excited 70 P3/2 atoms as 
a function of the pulse area for square laser pulses 
(iV = 70 and p = lO" cm"^). In (a) the pulse du- 
ration was varied and in (b) the single-atom Rabi 
frequency. In both graphs the dashed lines repre- 
sent some particular random distributions of atoms, 
while the solid lines are the average dependencies 
over 100 random arrangements of atoms. 



start the recursion. In the first step, we calcu- 
late the interactions between superatoms using 
Eq. (HI). In the second step, we check if the num- 
ber of superatoms is equal to the desired number 
(chosen in advance) of superatoms. If it is, the 
recursion is over, otherwise, we execute the third 
step. In this step, we group superatoms cor- 
responding to the largest |fcy|. Therefore, the 
number of superatoms is reduced by one and 
thus we need to recalculate all kij. This means 
that the recursion cycle is initiated again. 



Results and discussion 

We consider mesoscopic samples of '--^lO fj,m 
diameter containing up to 100 ground-state 
atoms. It is assumed that 5s ground state ru- 
bidium atoms are excited by one-photon transi- 
tions to 70p3/2 Rydberg states. The positions of 
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atoms within a sample are randomly generated 
before the time evolution of the system takes 
place. We typically consider samples of 70 atoms 
and density of 10^^ cm~'^. The only exception 
is the analysis of density fluctuations, shown in 
Fig [HI where the number of atoms varies accord- 
ing to a Poissonian distribution. The effect of 
changing the density is similar to changing the 
interaction strength, and so we can vary only 
one of them. We vary the scaled interaction 
strength which is obtained as follows. From 
the Schrodinger equation idip/dt = Hip, after 
the scaling i — > t/r, — > ilr and kij kijT, 
we conclude that the final excitation probability 
Poxc is a function of the pulse area and the prod- 
uct kijT ^ kij /T , i.e. Pcxc Pc-xci^'^ 7 kijT). 
Here r is the pulse duration (for a Gaussian 
it is the FWHM) and F is the bandwidth. In 
our simulations, the range of possible pulse ar- 
eas (~ fir) is the same for all types of laser 
pulses. Because the single-atom Rabi frequency 
n is scaled by in the regime where the 

correlations between atoms are significant, the 
many-body characteristics of Rydberg excita- 
tion should be present at relatively small fl. 

The effect of interactions can be amplified, 
not only by increasing kij, but also t. Increas- 
ing T also means reducing the bandwidth (for 
ideal pulses). When, for some r, the interac- 
tions between nearest neighbors become compa- 
rable with the bandwidth, the influence of inter- 
actions on the excitation dynamics should be- 
come noticeable. In our simulations the interac- 
tions between atoms are given by kij — Cq/ R^j . 
This type of interactions is also supported by 
the experiment p^ . in which the production 
of ions in collision processes is consistent with 
the assumption of an attractive van der Waals 
potential and a theoretical estimate of its mag- 
nitude. In the model presented in we cal- 
culated the effective vdW coefficient C% to be 
7C6/60, where Cq is the dispersion coefficient 
of the strongest np+np potential [ISj]. This Cq 
is essentially mo del- dependent, however, using 
this value we had a good agreement with the 
experimental data [4j] and thus we use it in this 
calculation as well. 

In Fig. [U the time and omega dependen- 
cies of the excitation probability are shown for 
square laser pulses in panels (a) and (b), respec- 
tively. These two dependencies are presented as 
functions of the pulse area. In both graphs the 
dashed curves correspond to particular random 



distributions of atoms within the sample. The 
solid curves are the average over the results ob- 
tained from one hundred random arrangement 
of atoms. If we use fifty different placements 
instead of one hundred, the difference in the av- 
erage excitation probability is only about 1 %. 
The evaluation of the omega dependence is much 
more demanding because for each value of f2 we 
have to calculate the corresponding time depen- 
dence first. Although excitation probabilities for 
particular random distributions of atoms clearly 
exhibit many-body oscillations, when averaged 
over many distributions, these oscillations are 
significantly suppressed. Collective oscillations 
presented in Fig. [3] are more pronounced be- 
cause a much greater scaled interaction strength 
is used. 

In Fig. [2] the average time and omega depen- 
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(a) Time dependence 
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Gaussian pulse 
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(b) Omega dependence 
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Pulse area (x Jt / N^'^) 

FIG. 2: The average dependencies of the number of 
excited 70 P3/2 atoms versus the pulse area for dif- 
ferent types of laser pulses (A'' = 70 and p = 10^^ 
cm~'^). In the absence of interactions all these de- 
pendencies should be the same. In (a) the pulse du- 
ration was varied and in (b) the single-atom Rabi 
frequency. Both solid lines correspond to square 
pulses and dashed lines to Gaussian ones. This plot 
demonstrates that the excitation blockade is a bit 
more efficient for Gaussian pulses. 
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dencies of the excitation probability for square 
and Gaussian pulses are presented. Here we 
show how different pulse shapes affect Rydberg 
excitation. It is known that for resonant single- 
atom excitation, the excitation probabilities de- 
pend only on the pulse area. However, from 
-fcxc = Pexc{^T,kijT), wc scc that varying the 
pulse area by changing r or SI leads, in general, 
to different results. Only in the limit of a full 
blockade, according to Eq. ([1]), does P^xc again 
become a function of the pulse area only. The 
average time dependencies, for both types of 
pulses, slowly decrease after reaching saturation 
(the first maximum). This slow decrease can be 
explained by the reduction of the bandwidth F 
for longer pulse durations, i.e. the blockading 
due to interactions becomes more effective. The 
slow increase with omega, after reaching satura- 
tion, is expected to happen due to higher laser 
power. This figure shows that the excitation 
probabilities for Gaussian pulses are systemat- 
ically lower than those for square pulses. This 
is also expected because of the following some- 
what oversimplified argument. From the point 
of view of mean-field theory, the effect of inter- 
actions is a level shift, which can be expressed as 
some effective detuning of the excited level. Ig- 
noring the time dependence of such level shifts, 
and in the first approximation, the probabilities 
are given by the Fourier transform of the pulse 
envelope. However, the Fourier transform of a 
Gaussian is also a Gaussian so that the Fourier 
components of a Gaussian vanish much faster for 
large detunings than the corresponding Fourier 
components of a square pulse (which fall off as 
~ 1/A). 

In the previous figures we have assumed the 
interactions between 70p3/2 Rydberg atoms and, 
for square pulses, pulse durations of 10 ns at 
most. However, the product kijT can still be 
significantly increased, either by going to higher 
principal quantum numbers n (increasing kij 
due to the n^^ scaling of Cg) or using longer 
pulses. In practice both ways should probably 
be used because for high n the diatomic energy 
levels become very close to each other and the 
interactions at short internuclear distances have 
complicated forms due to avoided crossings be- 
tween potential curves. However, increasing r 
also means decreasing F so that the excitation 
of atoms at shorter internuclear separations can 
be ignored and thus the van der Waals form of 
interaction can still be appropriate. Actually, 
the approximations in the superatom approach 



are justifiable if such excitation can be ignored. 
In fact, the largest Kpq have no infiuence on su- 
peratom interactions kij because for such Kpq, 
ideally, both atoms p and q should belong to 
the same superatom. The simulation shown in 
Fig. [3] is obtained for the product kijT 15 times 
larger than the value used in Figs. [T]151 This 
would roughly correspond to the interactions be- 
tween n = 90 Rydberg atoms with the same 
bandwidth F. Both curves are obtained after 
averaging over many random spatial distribu- 
tions of atoms. The oscillations are obviously 
more pronounced here than in Fig. [T][2l To see 
how robust these oscillations are against density 
fluctuations, we varied the number of atoms in 
the sample according to a Poissonian statistics. 
The dashed curves is obtained after this addi- 
tional averaging over density fluctuations. 

In the flrst approximation, we can find an ap- 
proximate parametrization of collective oscilla- 
tions as follows. From the fact that the initial 
excitation probabilities do not depend on inter- 
actions, and the assumption that as soon as the 
interactions start to dominate the process they 
quickly saturate the excitations, one readily de- 
rives the following approximate formula 

,x(lmax) p(lmax)2 

p(lmax) ^ ^''cxc ^ 

-^exc 4 ' 

where Po^^^^^ and Nci^^^^ are respectively the 
excitation probability and the number of excited 
atoms at the first maximum of Pcxc, and i^(i™ax) 
is the corresponding pulse area. The last for- 
mula is just the excitation probability for iso- 
lated (noninteracting) atoms if the pulse area is 
small. We can use it because, initially, there are 
no excited atoms and so there are no effects of in- 
teractions. In the next chapter we give a formal 
proof (??) for this claim. The effects of interac- 
tions for small pulse areas F are proportional to 
F'^. Depending on the strength of interactions, 
this saturation can be achieved much faster than 
the one in single-atom processes. In terms of 
frequency, the collective phenomena are faster 
than their counterparts in single-atom processes. 
Assuming that the phase (/fcoU of collective os- 
cillations can be simply characterized by some 
collective frequency recoil (i-e. (/icoii ~ f^coU'''), 
where f^coii = ctfl, the scaling parameter a can 
be obtained from the saturation condition 

^^(Imax) _ (g) 
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without fluctuation (N=70) 
with fluctuation (<N>=70) 



Qt (X It/N ) 



FIG. 3: The time dependence of the number of ex- 
cited atoms for the scaled interaction strengths kij 
15 times larger than those used in the other simula- 
tions (Figs. [l]and[2}. To check the robustness of the 
many-body oscillations against density fluctuations, 
the number of atoms in the sample is varied accord- 
ing to a Poissonian distribution. The solid curve is 
the dependence for Ai' = 70 and the dashed curve 
is the average dependence over density fluctuations 
with {N) = 70. 



Combining the last two equations one gets 




(7) 



To illustrate excitations in large samples with 
a significant excitation blockade, a domain (or 
"bubble") picture is often used. Each domain 
represents a region in which there exists ex- 
actly one Rydberg atom. Denoting the num- 
ber of atoms in a domain by Nd, then Nd — 



(linax) 



/N, a 



and 



(8) 



The factor ,/f = 0.886. 



1 can (and proba- 
bly should) be approximated by one (because for 
Nd = N, according to Eq. dJ), f^coii = VNfl). 
We conclude that it seems that \/Nd is the 
scaling factor of collective oscillations. We can 
verify this estimate by comparing with the nu- 
merical calculations shown in Fig. [2] and Fig. 
131 Also, if the phase of collective oscillations 
can be really described by a simple parameter 
^^coii (or y/Wo ) , then the ratio of the pulse ar- 
eas corresponding to the second and the first 
maximum of Pexc should be three. The figures 
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FIG. 4: The relative variance Ancxc/AnJ,xc of the 
number of excited Rydberg atoms versus the square 
of the pulse area for the same parameters used in 
Figs. IH2I The actual variance AtIoxc is expressed 
using the calculated variance Ariexo obtained from 
P, assuming no correlation between excited atoms. 
In (a) pulse duration was varied and in (b) the 
single-atom Rabi frequency. Both solid lines corre- 
spond to square pulses and dashed lines to Gaussian 
ones. 



clearly show that the actual many-body behav- 
ior is much more complicated, so we just want to 
verify how much it deviates from the simple pic- 
ture. For convenience, we define two parameters 

7 = a/y/NE and (3 = i^(2max)/^(lmax)_ jj^^gj._ 

estingly, for all square pulse results presented in 
Figs. and [31 we find /3 — 2.5, while for the 
Gaussian pulses in Fig. [2 /? = 2.3. The param- 
eter 7 varies more. For a square laser pulse in 
panels (a) and (b) in Fig. [21 7 is respectively 
0.95, 0.84 and 1.07. For the Gaussian ones in 
plots (a) and (b) in Fig. [1 7 is 0.87 and 0.77. 

Besides the many-body Rabi oscillations, 
there are other manifestations of correlations be- 
tween interacting excited atoms. One can study 
the fluctuations of the number of excited atoms. 
Without these correlations, the probability to 
have a certain number of excited atoms is deter- 
mined by the average excitation probability P 
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and the total number of atoms N. The proba- 
bihty P{ncxc) to get any number rioxc < of 
excited atoms is given by the BernouUi formula 



^'(nexc) 



N 



P"'=-(l-P) 



(9) 



For a given excitation probability P, one can 
calculate the expected variance Ang^^. assum- 
ing no correlations, and then compare it with 
the actual Arioxc- We take the relative size of 
these variances as the measure of these fluc- 
tuations. This procedure can be done experi- 
mentally as well. In Fig. [51 we show the ra- 
tio Anoxc/An^xj. as a function of the pulse area. 
This ratio is smaller than one because there are 
some restrictions on which combinations of ex- 
cited atoms can be likely excited in the presence 
of interactions. On the other hand, in the ab- 
sence of correlations, all combinations of excited 
atoms are equally probable. This figure also 
shows that after reaching the minimum, the ra- 
tio Anoxc/Arigxc increases again due to decoher- 
ence. The experimental results @ on counting 
statistics of somewhat larger samples than we 
consider here do show the sub-Poissonian char- 
acter of Rydberg excitation. 

We have also calculated the spatial correlation 
function between the central atom and other 
atoms in mesoscopic samples. For larger sam- 
ples, this function was calculated in [ll|. They 
also found some interesting correlations if there 
is some frequency chirp of the laser pulse. In 
the presence of chirp, there is a region of in- 
ternuclear separations R for which the correla- 
tion function is greater than one. In the absence 
of chirp this region does not exist. We do not 
include any chirp in our simulations because it 
is known that it destroys even the single-atom 
Rabi flopping. A delicate point for calculating 
the correlation function is that what we really 
calculate is the correlations beetwen superatoms 
and these superatoms are quite extended ob- 
jects. Their linear size is about one third of the 
sample radius. Averaging over many random 
arrangements of atoms improves the determina- 
tion of the correlation functions, but we do not 
insist on having many points for this function 
(Fig. [5]). In the supperatom approach, the cor- 
relation function c{p, q) between any atom p be- 
longing to superatom i and any atom q belong- 
ing to superatom j is the same for all p and q and 
equal to the correlation function C{i,j) between 
superatoms i and j. In addition, c(j),q) = 
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FIG. 5: The correlation function between a central 
atom and other atoms versus the distance between 
them for the same parameters used in Figs. IH2I 
This function is obtained by averaging over one hun- 
dred different random placements of 70 atoms. This 
means that 6900 atom pairs are used to get this de- 
pendence. 



for atoms p and q belonging to the same super- 
atom. This suggests that the correlation func- 
tion c{p, q) implicitly contains, in some sense, 
spatial averaging over the spatial extension of a 
superatom. Since, on average, excitation proba- 
bilities P{i) depend on the atom's location, the 
correlation function has to be 



c{p, q) 



P{p,q) 

P{p)P{q) ' 



(10) 



where P{p, q) is the probability to excite simul- 
taneously atoms p and q. In our case one of 
atoms p and q is always a central atom. The 
correlation function averaged over one hundred 
random positions of atoms is presented in Fig. O 
This figure shows that there are almost no cor- 
relations (i.e. c(p, q) ~ 1) between the central 
atom and the atoms near the sample surface. 



I. CONCLUSION 

We have investigated the many-body dynam- 
ics of Rydberg excitation in ultracold meso- 
scopic systems of ~ 10 ^m diameter. Vari- 
ous processes in Rydberg gases are studied in 
detail by numerically solving the many-body 
Schrodinger equation. The possibility of many- 
body Rabi oscillation of Rydberg excitation 
is explored assuming that the 5s ground-state 
atoms are excited to the 70p Rydberg state by 
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laser pulses of ~ 10 ns duration. For typical ex- 
perimental parameters, we have found that exci- 
tation probability for particular random distri- 
butions of atoms clearly exhibit many-body os- 
cillations. However, when averaged over many 
distributions, these oscillations are significantly 
suppressed. More robust collective oscillations 
could be achieved by going to higher Rydberg 
states n = 90, or decreasing the laser bandwidth 
(increasing pulse duration). We have evaluated 



the correlation function between a central atom 

and other atoms. It shows that excited atoms 
are strongly correlated within a range of few /xm. 
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